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The present work considers (4+l)-dimensional spatially homogeneous vacuum 
cosmological models. Exact solutions — some already existing in the literature, 
and others believed to be new — are exhibited. Some of them are the most general 
for the corresponding Lie group with which each homogeneous slice is endowed, 
and some others are quite general. The characterization "general" is given based 
on the counting of the essential constants, the line-element of each model must 
contain; indeed, this is the basic contribution of the work. We give two different 
ways of calculating the number of essential constants for the simply transitive 
spatially homogeneous (4+l)-dimensional models. The first uses the initial value 
theorem; the second uses, through Peano's theorem, the so-called time-dependent 
automorphism inducing diffeomorphisms. 
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1 Introduction 



Since the dawn of General Relativity people have been interested in finding exact solu- 
tions to Einstein's Field equations 1 . However, due to the fairly complicated nature of 
the field equations, we usually impose symmetries in order to make the field equations 
more tractable. Some of the most successful schemes of symmetry reductions are the 
so-called Bianchi models in (3+l)-dimensional cosmology [2-4]. Here, in this paper we 
will consider their (4+l)-dimensional counterparts [5,6]. 

The study of higher- dimensional models has - especially since the advent of String 
Theory [7, 8] - become increasingly popular in the recent years. For example, exact 
solutions like plane-wave spacetimes have been in focus for the last few years because 
they admit supersymmetry and they provide with an exactly soluble string background. 
Plane-wave spacetimes are some of the solutions of the models considered here. More 
specifically, we will consider (4+l)-dimensional spatially homogeneous spacetimes which 
are solutions of the vacuum field equations. Equivalently, we will consider Ricci-flat 
spacetimes; i.e. spacetimes obeying 

Rab = 0, (1.1) 

which admit a group acting simply transitively on the spatial hypersurfaces. The ques- 
tion we are addressing is "How large is the set of Ricci-flat spacetimes within the set of 
models considered?" , or equivalently "How many parameters are necessary, in principle, 
to specify a solution of this equation provided that they are spatially homogeneous?". 
The answer to this question turns out to be that 11 parameters are needed to be specified 
to give a solution for the most general classes. 

Many interesting phenomena are related to this issue. For example, a by-product of 
our analysis is that we are able to determine which are the most general vacuum models 
within the class of spatially homogeneous models. In (3+l)-dimensional cosmology, the 
most general simply connected Bianchi vacuum models, namely type VIII, IX, and the 
exceptional model VI* 1 ^ 9 [4] 2 , are all chaotic in the initial singular regime [9-18]. In 
4+1 cosmology one might wonder if the same is the case 3 . We will also give some 
exact solutions, thereby providing some examples of spacetimes of each class. In some 
cases, the entire family is known explicitly; in others only a few, or even none, are 
known. However, some of these special solutions have some interesting properties — like 
self-similarity — which may be important in the late-time behaviour for more general 
solutions (see e.g. [4]). As an explicit example of this, the plane- wave solutions - which 
will be discussed later - were shown in [20] to be the attractors within their class of 
models. 

The paper is organized as follows. Next, we introduce the automorphism group and 
see how it is related to coordinate transformations, or the gauge freedom in general 

1 See e.g. [1] for the known exact solutions in 3+1 dimensions. 

2 Strictly speaking, it depends how one counts. If we include the group parameter as an essential 
constant, then types Vlh and VII/j are equally general. 

3 This was recently addressed for some of the models in [19]. 
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relativity (see e.g. [21,22]). Then, in sectional we present our main results, namely 
the counting of the essential constants for the simply transitive, spatially homogeneous 
models of dimension 4+1. In sect ion 0] we present some exact solutions before we conclude 
in sectional 

2 The role of the Automorphism group 

Let us first exhibit some basic assumptions that lie in the foundation of our work. 

Spacetime is assumed to be the pair (J\4,g) where M. is a 5-dimensional, Hausdorff, 
connected, time-oriented and C°° manifold, and g is a (0, 2) tensor field, globally defined, 
C°°, non-degenerate and Lorentzian (i.e. it has signature (— , +, +, +, +)). In the spirit 
of 4+1 analysis, we foliate the entire spacetime like M. = MxS ( , where the 4-dimensional 
orientable submanifolds S t (surfaces of simultaneity) are space-like surfaces of constant 
time. The assumption of spatial homogeneity corresponds to imposing the action of a 
symmetry group of transformations G upon the manifolds S t . Usually, the group G is 
not only continuous, but also a Lie group — thus denoted by G r , where r is the dimension 
of the space of its parameters. Avoiding the details on these issues — these matters can 
be easily found in a standard reference, see e.g. [23, 24] — we simply state that spatially 
homogeneous models with a simply transitive action of the group G4 are described (apart 
from the topology of E t which we will assume is simply connected) by an invariant basis 
of one-forms u) a = erf (x)dx* (their Lie derivative with respect to the generators of the 
Lie group G 4 C a ( x )i are zero )- There is also the case of spatially homogeneous models 
in which the group acts multiply transitively (the number of generators is more than 4 
and there does not exist a proper invariant subgroup of dimension 4 acting transitively). 
The multiply transitive cases, which in dimension 4+1 are 5 in total [6], will not be 
considered here. More generally, a Lie Group, G r , is said to act transitively if the 
following requirements are satisfied: 

(i) r > dimension of 4) 

(ii) the rank of the matrix constructed by the generators (seen as vector fields) are 
everywhere equal to the dimension of S 4 (= 4). 

Geometrically, the above requirements imply that two different points in a given domain 
of S t can be interchanged by a Lie group transformation. Simply transitive action, which 
is our concern in this paper, corresponds to the case r = dimension of S t (= 4). In this 
case we note that the most general line element, manifestly invariant under the action 
of the group, takes the form (in an appropriate coordinate system): 4 

ds 2 = (N a (t)N a (t) - N 2 (t))dt 2 + 2N a (t)of (x)dtdx i + 7 ^(t)af (z)af(x)dxW, (2.1) 

Greek indices label the invariant one- forms while the Latin indices run over the spatial coordinates 
(i.e. from 1-4). 
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with 



°tj{*) ~ °U X ) = 2C >i ( 2 - 2 ) 

where "f a p(t) is the metric induced on the surfaces E t (and thus constant on them); N(t) 
is the lapse function; N a (t) is the shift vector (N a (t) = j a ^(t)N p (t), ^ al3 (t) being the 
inverse of 7 Q/ g (£)); and are the structure constants of the corresponding (closed) Lie 
algebra. In 4 dimensions, there are 30 closed, real, Lie algebras [25,26]. 

At this point, a question arises; is there any particular class of General Coordinate 
Transformations (G.C.T.s) which can serve to simplify the form of the line element and 
thus also Einstein's Field Equations (E.F.E.s)? The answer is positive and a thorough 
investigation of this problem and its consequences is given in [22]; indeed, not only is 
there a class of G.C.T.s which preserves the manifest spatial homogeneity of the line 
element (|2.1j) . but it also forms a continuous (and virtually, Lie) group. This group is 
closely related to the symmetries of the symmetry Lie group G>; it is its automorphism 
group. 

In the spirit of the 4+1 analysis we consider, apart from the time reparameterization, 
the following G.C.T.'s : 

t ^t = t^t = t 
x i _> & = g\t,x j ) <=> x l = f{t,x j ). (2.3) 

After insertion of ()2.3|) into ([2.1)1 . the wish to preserve the manifest homogeneity of the 
latter leads, in a first step, to the allocations: 

= <i<J)P a (t,x) (2.4a) 
^ = ^(f)A«(t,x)^(x), (2.4b) 

and, consequentlythe definitions 

N(t) = N(t) (2.5a) 
N a (t) = S^(t)(N p (t) + P p (t, x)) (2.5b) 
j a p(t) = A£(t, x)A%(t, x) lllv {t), (2.5c) 

where N a (t,x) = ^ a P{t)Np{t,x) with P l a {x) being the inverses of crf{x) — quantities 
which exist in the simply transitive cases. In order for the transformations ()2.3|) to have 
a well defined non-trivial action, it is pertinent for the quantities A^ and P a to be space 
independent. So, 

%■ = <U)P a {t) (2.6a) 
% = <U)m)°l^ (2-6b) 
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and therefore 



N(t) = N(t) (2.7a) 
N a {t) = S^(t)(N^(t) + P p {t)) (2.7b) 
7<*(t) = A»(t)A%(t)^ v (t). (2.7c) 

Thus (I2.f)j) instead of being allocations, turn into a set of highly non-linear partial 
differential equations. Integrability conditions for this system, i.e. Frobenious' Theorem, 
results in the system (the dot, whenever used, denotes differentiation with respect to 
time): 

C$J${t) = C&AfflAift) (2.8a) 
h^t) = C^(t)A^t), (2.8b) 

and "Time-Dependent Automorphism Inducing Diffeomorphisms" (A.I.D.s) emerge. The 
automorphisms of a Lie group G r form a continuous group. Those members of the group 
which are continuously connected to the identity element, form a Lie group as well 
- even though the topology of the latter might be different from that of the former. 
If one considers parametric families of the automorphic matrices, characterized by the 
parameters t\ Ag(t; t 1 ), and defines: 

AJgfor*) . =5<l (2.9a) 

t'=0 

dA a (t]T* 



dr { 



n = x h) ( 2 - 9b ) 



where are the generators with respect to the parameter r % of the Lie algebra of the 
automorphism group, then from the first of (|2.8j) . after a differentiation with respect to 
t\ one gets 

For an extensive treatment on these issues see [27], while for the relation and usage of 
these generators with conditional symmetries, see [28,29]. 

In the n + 1 decomposition of the spacetime (here n = 4), the E.F.E.s in vacuum 
assume the form: 

E° = K«KP-K 2 + R = (2.11a) 
K = ^C^-K^ = (2.11b) 
E% = K$ - NKKp + NB% + 2N'(KZC$ p - K»C a up ) = (2.11c) 
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with 



Kp{t) = r p (t)K p0 (t) (2.12a) 
B$(t) = l ap (t)R P p(t) (2.12b) 
Kafi(t) = -^L-fi^t) + 2 lau (t)C^(t) + 2 lf}v {t)Cl p N0{t)) (2.12c) 



+ 2C x aK C^ x {t)^{t) + 2C^C; iy7 aA(*)7^(0- 



(2.12d) 



Since G.C.T.s are covariances of the E.F.E.s, the same form of equations ()2.11|) holds 
for the transformed quantities: 

E° = E° = Q (2.13a) 

E% = AZE% = (2.13b) 

^ = S£Aj££ = 0, (2.13c) 

where is the inverse of A^.This can be explicitly seen by observing that the extrinsic 
curvature transforms as a (0,2) tensor under these transformations, despite the mixing 
of time and space coordinates. The effect of a time reparameterization is trivially seen 
also to be a covariance . Finally some terminology is needed; ()2.11a|) is called "Quadratic 
Constraint", ()2.11b|) are called "Linear Constraints", and (|2.11c|) are simply the "Equa- 
tions of Motion". 



3 Essential Constants 

The task of finding the maximal number of essential constants for each model is compli- 
cated by the presence of the quadratic and linear constraint equations. 
The first thing to observe is that their time-derivatives vanish by virtue of the spatial 
equations of motion; therefore, they are first integrals of motion for these equations and 
they will be satisfied at all times once they are satisfied at one instant of time. The 
constraint equations can thus be considered as algebraic relations restricting the initial 
data at some arbitrarily chosen hypersurface. Accordingly, one initial datum will be 
absent for each such independent constraint. 

The second important thing is that the additive constant of integration at the right-hand 
side of these constraint equations is identically zero. This points to the fact that the 
presence of these equations signals the existence of "gauge" symmetry for the whole sys- 
tem of equations, namely the time-dependent A.I.D.'s briefly described in the previous 
section. Under these transformations, one more constant becomes absorbable. Thus, if 
we wish to consider the constraints as full fledged (first class) differential equations, we 
have to subtract two degrees of freedom ( constants in our case) for each such indepen- 
dent equation. 

Both points of view are correct and valid: they are nothing but different aspects of the 
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same ingredients of the theory of differential equations. Thus they should yield the same 
final result concerning the maximal number of essential constants. Bellow we present 
the counting algorithms of this number for all 30 (4+1) simply transitive, spatially ho- 
mogeneous vacuum geometries, according to both points of view. 

3.1 The Initial Value Theorem 

In this section, we apply the initial value theorem to find the maximal number of essen- 
tial constants each line element should contain in order to describe the entire space of 
solutions for the given model. In 3+1 dimensions, such a counting has been done some 
time ago (see e.g. [4] and the older references therein) using the Behr decomposition of 
the structure constants C^ 7 for 3-dimensional Lie algebras. However, such a decomposi- 
tion is not known for Lie algebras of dimension 4 or higher. Hence, we have to apply an 
alternative counting procedure in order to find the essential constants. 

A counting which is independent of the dimension can be given using the initial value 
theorem. This theorem is stated for the (3+l)-dimensional case in, for example, Wald's 
book [30]. However, it is fairly easily seen that this theorem is valid in any dimension; 
the arguments in the proof does not depend explicitly on the dimension of the spacetime. 

Roughly speaking, the initial value formulation says that a spacetime satisfying the 
Einstein equations is uniquely determined by specifying the metric, h a p, and the corre- 
sponding extrinsic curvature, K a @, of an initial spatial hypersurface (i.e. ja^to) = h a p) 
-at the Gauss normal coordinates system, in which the shift vanishes. The initial 
data must also satisfy the quadratic constraint, and the linear constraint on the ini- 
tial hypersurface which are purely algebraic in the initial data. Furthermore, isomet- 
ric diffeomorphisms on the initial hypersurface, can always be extended to isometric 
diffeomorphisms of the entire spacetime. 

The theorem does not mention whether two different initial data can lead to the 
same spacetime. However, any initial data always generates a one-parameter family of 
data which will yield the same maximal development. This one-parameter family is 
exactly the time evolution of the pair ( r y a p(t), K a p{t)). Hence, for a spacetime foliated 
into spatial hypersurfaces, any hypersurface may serve as an initial hypersurface. 

The initial value formulation thus provides us with the following algorithm for count- 
ing the essential constants for the spatially homogeneous model of type A: 

#(h a p, K a p) — dimAut(A) — ^(independent constraints) — 1. (3.1) 

In our case, #(h a p, K a p) = 20, since h a p and K a/3 are symmetric 2-tensors. Ant(A) is 
the automorphism group for the Lie algebra A; these automorphisms can be seen to be 
the effect of isomorphic diffeomorphisms on the initial hypersurface. Thus they carry the 
relevant "gauge" freedom which must be subtracted [27,31]. On the initial hypersurface 
the constraints (quadratic plus linear constraints) are only algebraic equations, thus 
subtract one for each constraint. Finally, we subtract 1 due to the fact that each initial 
hypersurface traces out a one-parameter family of initial data each giving rise to the 
same spacetime. 
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Using the above algorithm we produced Tables 121 and El giving the number of essen- 
tial constants for all 30, simply transitive, spatially homogeneous vacuum cosmological 
models of dimension 4+1. Table El contains the essential constants for the general form 
of the algebras, while Table 01 contains the essential constants for the exceptional cases in 
which for some values of the parameters (of the Lie algebra) some of the linear constraints 
vanish identically. 

3.2 Time Dependent A.I.D.s 

In this section, we apply the time dependent A.I.D.s to perform a second independent 
counting of the maximal number of essential constants each line element should contain 
in order to describe the entire space of solutions for the given model. This way of counting 
is valid in any spatial gauge. The key observation is that the solutions to the integrability 
conditions (J2.8)) always contain 4 arbitrary functions of time. These arbitrary functions 
are distributed in and P a in a way that differs for each of the 30 models; e.g. to 
take an extreme case in the Kasner-like model, 4A\, Ajg is completely constant while all 
4 arbitrary functions of time are located in P a . In all cases, P a contains all arbitrary 
functions through either their derivatives, or themselves. Thus two distinct ways of using 
the gauge freedom suggest themselves, leading to two versions of the counting algorithm: 
The first, is to use the whole freedom in order to set the shift N a equal to zero and 
then see how many first class linear constraints remain. The corresponding version of 
the algorithm is: 
D = 2x (# of 7a/3 ) 
—2 x # ( linear constraints ) 
— 2x (the Quadratic Constraint) 

(parameters of Outer Automorphic matrices) 

-(«) 

where k = dim(Inner) — # functionally independent Linear Constraints. 
Peano's theorem requires 2 initial data for each 7 Q/ 3 since the system is of second order. 
We subtract 2 constants for each independent first class constraint. Finally we subtract 
the remaining rigid symmetries which are the parameters of the outer Automorphisms 
plus the difference between the number of parameters of the inner automorphisms sub- 
group and the number of functionally independent linear constraints. 

The second consists of all other options, e.g. we can use the functions of time con- 
tained in to simplify the scale factor matrix 7 Ql a(£) and the remaining functions 
contained in P a -if any- to alter somehow the initial shift vector (e.g. equating some 
components or setting some of them equal to zero). Now the algorithm reads: 
D = 2 x of 7 a/ g) + 1 x of possibly remaining 5 shift vector's components) 
—2 x # (of those linear constraints which do not finally involve shift vector's compo- 
nents) 

5 i.e. after solving algebraically as many as linear constraints is possible -in terms of the shift's vector 
components-, i.e. the shift components which are not expressed in terms of the scale factor matrix 
components and their derivatives. 
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— 2x (the Quadratic Constraint) 

— # (parameters of those Outer Automorphic matrices which preserve the form of the 
reduced 7 aj g) 

For the sake of illustration we give bellow three examples of counting with both 
versions of the algorithm presented in this subsection. 



3.2.1 Type A 2 © A x 

The structure constants are: Cf 2 



A?(f) 

Kit) 



1. Thus: 
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1 











\ 




A 5 (0 


A 6 (0 












A 9 





An 


Al2 




V 


Al3 





Al5 


Al6 


/ 



P a {t) = 



(\ 6 (t) \' 5 (t)-Ht) m) 



2A 6 (t)' 2A 6 (f) 

Four functions of time appear -as expected; two of them in A^(t) and correspond to the 
inner automorphism proper invariant subgroup. The number of functionally independent 
linear constraints, is 4. 



1st version We use our entire freedom in order to set the shift vector equal to zero. 
So: 

^ Linear Constraints in terms of 7^ = 4 
# of parameters of Out. Aut. matrices = 6 

k = 2 - 4 = -2 
D = 2x 10 - 2x4-2-6 - (-2) = 6 



2nd version We use our freedom in order to set: iV 3 (t) = N 4 (t) = and 7 12 (t) = 0, 
722 (t) = 1. So: 

# 7a/3 = 8 

# remaining iV a = 
# Linear Constraints in terms of 7q/ 3 = 2 
# of parameters of those Out. Aut. matrices which preserve the form of 7 Q/ g = 4 

5=2x8-2x2-2-4=6 
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3.2.2 Type A 3i6 © A l 

The structure constants are: Cf 3 = — 1 and C23 



1. Thus: 



AS(f) 



/ ccos(/(t)) csin(/(0) A 3 (t) \ 

-csin(/(t)) ccos(/(t)) A 7 (t) 

10 

V A15 Ai 6 y 



P a (t) 



(a 3 (o /'(o) - a' 7 (o -(\ 7 (t)f(t))+m -f\t) 



Four functions of time appear -as expected; the three in A^(t) and correspond to the 
inner automorphism proper invariant subgroup. The number of functionally independent 
linear constraints, is 3. 

1st version We use our entire freedom in order to set the shift vector equal to zero. 
So: 

# la(3 = 10 

# Linear Constraints in terms of 7 Q/ 3 = 3 
# of parameters of those Out. Aut. matrices which preserve the form of 7 Q/ g = 3 

K = 3-3=0 
D = 2x 10 - 2x3-2-3 = 9 



2nd version We use our freedom in order to set: N 1 ^) = N 2 (t) = N 4 (t) = and 
712(f) = 0. So: 



# 7a/3 

# remaining N a 
# Linear Constraints in terms of 7^/3 
# of parameters of those Out. Aut. matrices which preserve the form of 7^/3 

£) = 2x9-2x2-2-3 



9 

2 
3 
9 



1 _i 

3' 3 



3.2.3 Type A 4 5 

The structure constants are: C\ A = 1, C| 4 = — | and C| 4 = — |. Thus: 



AS(0 



/ Ax(*) 



V 





el 
Ai(t) 3 

c3 
Ai(t)S 







c2 
Ai(t)3 

c4 

Al(t)3 





A 4 (t) \ 
A 8 (t) 
A 12 (0 
1 / 
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pn ( - (\,(t) K(t) - X^t) \' 4 (t)) ^ -(A 8 (t) AU^ + SAx^A^t)) ^ g ^ -Ai(*) 



2Ai(t) ' 2Ai(t) ' ' 2Ai(t) J 

Four functions of time appear -as expected; all the four in Ajg(t) and correspond to the 
inner automorphism proper invariant subgroup. The number of functionally independent 
linear constraints, is 2. 

1st version We use our entire freedom in order to set the shift vector equal to zero. 
So: 

# la(3 = 10 

# Linear Constraints in terms of 7 Q/ g = 2 
# of parameters of those Out. Aut. matrices which preserve the form of 7 Q/ g = 4 

K = 4-2 = 2 
D = 2x 10 - 2x2-2-4-2 = 8 



2nd version We use our freedom in order to set: 7n(t) = 1 and 714 (t) = 724 (t) = 
734 (0 = 0. So: 

# 7a/3 = 6 

# remaining N a = 2 
# Linear Constraints in terms of 7 a/ 3 = 
# of parameters of those Out. Aut. matrices which preserve the form of 7^ = 4 

£1 = 2x6 + 2- 2-4 = 8 



4 Exact solutions 

We will here provide with examples of spatially homogeneous vacuum solutions in 4+1 
dimensions 6 . There are some general things worth noting. For the decomposable cases, 
A 3 © Ai, we can generate vacuum solutions from scalar field solutions of the Bianchi 
models in 3+1 dimensions. More explicitly, given a vacuum solution in 4+1 dimensions 
with metric 

ds 2 5 = ds\ + e^dy 2 , (4.1) 

the metric ds\ = e~^ds\ will be a solution to the (3+l)-dimensional Einstein equations 
with a scalar field. Thus, by going the other way, we can construct vacuum solutions in 
4+1 dimensions from scalar field solutions in one dimension lower. In many cases (like 
type VIII©R and IX©R) these are the only non-trivial solutions one knows explicitly 
(see [!])■ 

6 Some of the solutions are previously known, even though in many cases the true number of free 
parameters was not recognized. 
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The main object of this section is to give some examples of solutions of the various 
types. For only two types we know all the possible exact vacuum solutions, the remaining 
cases we only know some special ones. 



4.1 4A 1 =I©R 

There is a 2-parameter family of Kasner solutions which exhaust all solutions of this 
type [32,33]: 

4 

ds 2 = -dt 2 + * 2pi dxMx\ (4.2) 

i=l 

where ^Pi = EiPi = 1 - 

4.2 A 2 2Ai =III©R 

There is a 2-parameter family of plane- wave solutions which can be obtained by restrict- 
ing the type VI^ © R plane-waves (III=VI_i) (see section l4~£|) . 

Also, there is a 2-parameter family of solutions, with a higher symmetry, given by: 

p 2 -2(l+a)i j f 2 p -2(l+o)t 

ds 2 = - 6 , — + ^— 2 (dx 2 + e-^dy 2 )+e 2at dz 2 + e 24 dw 2 ,(4.3) 

sinh cut sinh oot 

u 2 = a 2 + a + l. 

The symmetry group of these solutions is G5 = SL(2, R) x R 2 which acts transitively on 
the spatial hypersurfaces. 

4.3 2A 2 

There is one solution which can be obtained by a Wick rotation of a solution in [35] : 

ds 2 = -dt 2 + - [(dx 2 + e 2l 'dy 2 ) + (dz 2 + e 2z dw 2 )] . (4.4) 
3 

This has indeed the bigger symmetry group G$ = SX(2,R) x SX(2,R) acting on the 
spatial hypersurfaces S t . It is also algebraically special of type 22 in the sense of [35]. 

4.4 A 3) i Ai =II0R 

The general solutions (containing 6 parameters) are not known to our knowledge, but 
we have found a 4-parameter family of solutions. 7 It is given by: 

rfs 2 = _^ e (a 1+ a 2 +3a 3 )t cosh ^ dt 2 + ^_!^!L (dx _ zdy) 2 

ijj 04 cosh cut 

+e a 3 t^^l ( e -i*dy 2 + e a2 *dz 2 ) + e 2a3 'dw 2 , (4.5) 
uj v ' 



r See also [34] which considers the ^3,1 A\ and A 3j3 © A\ cases. 
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where u 2 = a^a 2 + 2(ai + a 2 )a 3 + a|. 

This family of solutions generalizes Taub's type II vacuum solutions. 



4.5 ^3,2 e Ai =iveR 

A 3-parameter family of plane- wave solutions is given by eq. (J4.7|) with s = 2j3 + . 

4.6 ^3,3 e Ai =veR 

A 2-parameter family of plane- wave solutions is given by eq. fl4.15j) with s = 2f3+. 
There is also 3-parameter family of solutions given by: 

ds 2 = 5 h e 2ait dw 2 

4 sinh d tut 

+ lfe^ ( e a 2 t e -2z dx 2 + e -a 2 t e -2z dy 2 + ^ ^ 

3u 2 = ?>a\ + a 2 . 

4.7 A 3A Ai =VI e R 

There is a 1-parameter family of solutions given by eq. 1)4.19)1 with p = —1, and q = 0. 

4.8 A\ h e Ai =vi A e r 

A 3-parameter family of plane- wave solution is given by eq. (|4.15j) with s = 2(3+. There 
are also a 2-parameter family of solutions given by eq. ()4.18|) with q = 0. 

4.9 ^3,6 e Ai =vn e r 

Apart from the solutions with higher symmetry deducible from scalar field Bianchi type 
VIIo, the authors do not know of any other non-trivial solutions. 

4.10 Af )7 © Ai =VII^0R 

A 3-parameter family of plane- wave solution is given in eq. (|4.2fij) with s = 2j3 + . 

4.11 A 3j8 Ai =VIII0R 

Apart from the solutions with higher symmetry deducible from scalar field Bianchi type 
VIII, the authors do not know of any other non-trivial solutions (see also [36]). 
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4.12 ,4 3 ,9 e ,4i =IX0M 

Apart from the solutions with higher symmetry deducible from scalar field Bianchi type 
IX, the authors do not know of any other non-trivial solutions (see also [36]). 

4.13 A 4 ,i 

No vacuum solutions of this type is known to the authors. 8 



4.14 A\ 2 

There is a 3-parameter family of plane- wave solutions for p > —2 [38]: 
ds 2 = 



,2t ( _ dt 2 + dw 2 )+e 2s( W+ t) 



X 



+ e 



20 + (w+t) 



(dy + Q 3 (w + t)dz) 2 + e 2Mw+t) dz 2 



(4.7) 



where 



and 



s(l-s) = 2P 2 + + l(Ql + Ql + Qt) 
Pi = 3/3+, 

3(3 + Q 2 - QiQ 3 



A = 



3p. 



+ 



R Q1Q3 



(4.8) 



(4.9) 



The group parameter is given by: 



V = 



(4.10) 



For p = —2, there is a 1-parameter family of vacuum solutions due to Demaret and 
Hanquin 9 [37]: 



ds 2 = k 2 e 3t2 t~^ (H + (-dt 2 + dw 2 ) + tle^dx 2 

_ i 

+tS(H + t^)e- 2w dy 2 + _\ 3 i e - 2u, (dz - u>dy) 2 . 

t 2 + t2 



(4.11) 



8 There are some known self-similar solutions with a perfect fluid [6]. Note there is a typo in eq. (83); 
all exponents should be divided by 7. 

9 However, they did not realize that the solution was part of a one-parameter family of solutions. 
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4.15 A\ 2 

There is a 2-parameter family of plane-wave solutions of one sets Q\ 
(3+ = in eq. lj£7|) . 



0, and then 



4.16 A 



4,3 



Plane-wave solutions for the Lie algebra type can be obtained by taking the p — > oo 
limit of ^4 2- in this limit we get f3 + = —s and thus the metric can be written [38]: 

ds 2 = e 2 '(-dt 2 + dw 2 ) 

+ e 6s(u,+ * ) (d* + ^i e - 3s ^+*)dy + [A + £(w + t)} e- 3s{w+t) dz^j 

+ (dy + g 3 (w + t)dz) 2 + dz 2 , (4.12) 



where 

* = I (i ± y / i-2(Q 2 + g^ + Q!)) , 

and A, 5 are given in eq. (|4.9|) with /3 + = — s. 
4.17 ^4,4 

There is a 3-parameter set of plane- wave solutions given by [38]: 

ds 2= e 2t ( _ dt 2 + dw 2 )+e2s( ^) 



(4.13) 



x 



dx + Qi(u; + t)dy + (w + t) 



Q2 + 



(w + t) 



dz 



where 



4.18 A% 



+ (dy + Q 3 (u; + t)dz) 2 + dz 2 



1 



(4.14) 



S (1- S ) = -(Q 2 + Q 2 + Q 2 ). 



Given p + g+ l>0,a 3-parameter set of plane- wave solutions can be given by [38]: 

ds 2= e 2 t( _ dt 2 + dw 2 ) + e 2 S ( W+t ) 

-4/3+ (uH 



V Pi P3P2 J 

+ e 2(f3 + +V3P~)(w+t) ^ dy + 03 e P-4w+t) dz ^ + e 2(J3 + -V3l3-)( W +t) dz 2 



(4.15) 
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where 



s (l- s ) = 2(/3 2 +/3 2 _) + |(Q 2 + Q 2 + Q 2 ) 
Pi = 3/?+ + v/3/3_ , P 2 = 3(3 + - v/3/3_ , P 3 = -2 V3/3- 

The group parameters are related to these parameters as follows: 

_ s + (P+ + V3P-) _ s-2/3 + 
P ~ s + {f3 + - q ~ s + (P + - y/3P_)' 



(4.16) 



(4.17) 



There are also some other solutions due to Demaret and Hanquin [37] 10 . Given 
p + q + 1 7^ 0, then there is a 2-parameter family of vacuum solutions: 



ds 2 = fc 2 (sinht) 2 £ p * (tanh0 2EQlPl (-dt 2 + dw 2 ) 



+ 



^(sinht) 2P * ftanh-Y^e 2 ^ (dx<) 



(4.18) 



where £P< = 1, £oj = 0, and £o: 2 = 1 + E^f- 

Given p + g + 1 = 0, there is a 1-parameter family of solutions due to Demaret and 
Hanquin [37]: 

ds 2 = fcV 1+p2+,?2 ^H(-dt 2 + dw 2 ) + ti (e 2w dx 2 + e 2pw dy 2 + e^dz 2 ) .(4.19) 

Also, for the exceptional case (q = — (1 + p)/2), we have found a self-similar 
solution which generalizes the Collinson- French type VI* wg vacuum: 



ds 2 



-dt z + £ 2 dx 2 + 



£ ( b } exp (— \/6(l + p)rx") dy H — tdx 

^ 7 * 2r\/b 



6(i+p) / i— \ 2 6p(i+p) / y— \ 2 
+t 6 exp ( 4v6rx J dz + £ 6 exp ( 4pv Qrx J dw" 



where 



a/1 + p + p 2 
5p 2 + 2p + 5 ' 



b = 5p 2 + 2p + 5. 



(4.20) 



(4.21) 



4.19 Aft 

This is a special case of the above. For the plane- wave solutions, eq. (|4.15|h one has to 
set Q\ = 0, and then P\ = 0. 



10 



They only give it as a 1-parameter family. 
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4.20 A% 

Similarly as in the above case, but now set Q 2 = and then P 2 = in eq. (|4.15|) . 

In addition to this, we have found a 3-parameter family of solutions for the particular 
value p — — 1 . It is given by: 



ds 2 



8u z 



a; 2 fc 2 e -(4a 1 +2a 2 )t dt 2 

- + 



2u: 



-a2*J_2 



sinh uit 



sinh uit 



dx" + dz 2 



L2 -(4ai+2a 2 )tJ w 2 
+e (2a 1+ a 2)ie 2» ^2 My 2 + *J dW 

sinh cut 

3a 2 + 3aia 2 + a 2 ,. 



(4.22) 



4.21 A\j 

The whole set of solutions is in this case known 11 . The set is 2-dimensional and the 
general solution is given by eq. ()4.18|) with the restriction Pi = P 2 = P3 = 1/3. 
Explicitly, 



ds 2 = k 2 sinhn (-dt 2 + dw 2 ) 

+ smh^te^ w (tanh- 
V 2 



2(4! 



dx + ( tanh 



t 



2a 2 



dy 2 



3a 2 + 3a 2 , + 3aia 2 . 



( tanh 
V 2 



, \ -2(ai+a 2 ) 



dz 2 

(4.23) 



4.22 Af fi 

Again we have plane- wave solutions [38]: Let s be given by 

s(l -s)= 2f3 2 + + ^lu 2 sinh 2 2(3 + \(Q\ + Q\). 



6 



Define also the two one-forms: 



uj 2 = cos[lj(w + t)]dy — sm[uj(w + t)]dz 
uj 3 = sin[ui(w + t)]dy + cos[ui(w + t)]dz. 

The plane-wave solutions of type A^ Q 6 can now be written: 

ds 2= e 2*(_ dt 2 + dw 2 )+e 2^ +i ) 

x e -^ w+ ^ {dx + e 3 ^ w+t ^ ( gi e-^ 3 - q 2 e f3 u 2 ) }' 
+ e 2 ^ w+t Ue- 2l3 (u 2 ) 2 + e 2f3 {u: 3 f 



L All solutions with a 7-law non-tilted perfect fluid is also known, see [6]. 



(4.24) 



(4.25) 



(4.26) 
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where 



qi= u* + 90l ' Q2= ^ + 9ft • (4 ' 27) 



The group parameters are related to these constants via 



4.23 A 4>7 

No such solutions are known to the authors. 

4.24 A Aj8 

There is a solution which is the p — > — 1 limit of the metric ()4.29j) . 



4.25 A% 

There is a simple power-law solution for each — 1 < p < 1: 



ds 2 = -dt 2 + t 2 dw 2 + k 2 t 3(p a +p+D e - 2(p+1)< ™(dx - zdy) 2 

2(p+2) 2 2(2p+l) 2 
+t 3( P 2 +P+ i) e -2^ dy 2 + t3(p 2 +p+1)e -2p<T«, dz 2 ; 



where 



2 = 7p 2 + 13p + 7 2 = 2(7p 2 + 13p + 7) 
°" 6(p 2 +p + l) 2 ' 9(p 2 +p+l) ' 1 } 

Due to the power-law dependence, this solution is self-similar. 

There is one special case worth noting, namely p = —1/2. 12 In this case the metric 
can be written: 

ds 2 = -dt 2 + ^ 2 (dw 2 + e- 2w dy 2 )+e- u '(dx-2dy) 2 + e u 'dz 2 . (4.31) 

Note that the spatial hypersurfaces are fiberbundles over H 2 . In fact, the symmetry 
group is larger for this metric than one would expect; it is the semi-direct product 
G 5 = R 2 x SL(2,R) with a U{1) stabilizer. 

12 This corresponds to a Lie algebra acting simply transitive on the model geometry F 4 [39]. 
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4.26 A{ 9 

There is a solution obtained from eq. ()4.29j) by setting p = 1, which is fairly interesting 
]. By a rescaling of the coordinates the solution can be written: 



ds 2 



+ 



dw 2 + e~ 2w 



dx + ^(ydz-2dy) ) fe^dy 2 + dz 2> 



(4.32) 



In this case the spatial surfaces are isometric to the complex hyperbolic space, H^, and 
hence, it has an 8-dimensional isometry group, Gs = PU(2, 1), acting multiply transitive 
on the spatial surfaces (it has a U{2) stabilizer). 



4.27 A% 

There is a solution obtained from eq. ()4.29j) by setting p = 0. 



4.28 A 4 ,io 

This algebra acts simply transitive on Nil 3 x R [39], so all solutions of the type II©K. 
admitting an extra symmetry acting on the spatial surfaces, are also invariant under 
this algebra. Hence, the solutions (|4.5jl with a\ = are invariant under this group. 
Solutions with A 410 as a maximal symmetry are not known to the authors. 



4.29 A\ 1X 

The solution ()4.32|) is invariant under this group due to the fact that this algebra acts 
simply transitive on [39]. Other solutions are not known to the authors. 



4.30 A Ajl2 

This algebra acts simply transitive on I 3 x R [39] so all solutions having this higher 
symmetry group are invariant under A 412 . An interesting example - although far 
from general - is the 1-parameter family of solutions obtained by Wick-rotating the 
5D Schwarzschild solution: 

ds 2 = + i (t 2 + 2M) dx 2 + t 2 [dy 2 + e- 2 ^(dz 2 + dw 2 )] . (4.33) 

As is clearly seen, this solution has a far larger symmetry group than A^\2, namely 
G-j = SL(2, C) x R. However, solutions with a maximal symmetry group A^u are not 
known to the authors. 
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TABLE 



Lie Algebra 


Non Vanishing Structure Constants 


4Ai 




A 2 © A x 


C? 2 = 1 


2A 2 


C 12 = 1 C 34 = 1 


A 3 ,i © A x 


C23 = 1 


A 3j2 © Ax 


C 13 = 1 C 23 = 1 C 23 = 1 


A 3 ,3 © Ax 


^13 = 1 ^23 = 1 


A 3 , 4 © Ai 


C 13 = 1 C 23 = — 1 


A% 5 © Ai < 


a < 1 


Cl3 = 1 ^3 = Ot 


A 3 , 6 © Ai 


c\ 3 = -1 c] 3 = 1 


A£ 7 © Ax < a 


C\ 3 = a Cf 3 = -1 C\ z = 1 C| 3 = a 


A 3 , 8 © Ax 


C23 = 1 ^13 = — 1 C?2 = — 1 


A 3 , 9 © Ai 


G 12 = 1 G 23 = 1 G 31 = 1 


A 4 ,i 


o 24 = 1 o 34 = 1 


^4,2 a ^ (0, 1) 


C 14 = a C A U = 1 C| 4 = 1CJj = 1 


^1,2 


G 14 - 1 G 24 - 1 G 34 - 1 G 34 - 1 


A 4 ,3 


/nrl 1 /"~(2 1 
U 14 — 1 L/ 34 — 1 


A 4 ,4 


o 14 — 1 o 24 — 1 o 24 — 1 o 34 — 1 o 34 — 1 


A 4 f -1 < a < < 1, a/3 ^ 


r~<l 1 /nr2 r~i?> a 
o 14 — 1 o 24 — a o 34 — p 


A 4 ; 5 ° -1 < a < 1, a ^ 


ri'Y 1 /~i2 „, /~i'S „, 

w-^ 4 — J. ^- / 24 — 3 4 — 


A^ 1 -1< a < 1, a ^ 


A 1 l -1 / * 2 /~~f3 1 

(_/ 14 — 1 o 24 — a (_/ 34 — 1 


^4,5 


14 = 1 24 = 1 L/ 34 = 1 


A^ a ^ 0, P > 


C! 4 = a C| 4 = /3 C 2 3 4 = -1 C| 4 = 1 C| 4 = /3 


A 4J 


G 14 = z G 24 = 1 G 34 = 1 G 34 = 1 G 23 = 1 


A 4 , 8 


Q23 = 1 c| 4 = 1 C34 = — 1 


<9 < 1/9 


< 1 


Q>3 = 1 Ci4 = 1 + /3 C| 4 = 1 C| 4 = f3 


^4,9 


^23 = 1 ^14 = 2 C| 4 = 1 C| 4 = 1 


^4,9 


C*23 = 1 C*14 = 1 ^24 = 1 


^4,10 


^23 = 1 ^24 = — 1 C34 = 1 


A° n a > 


G^ = 1 G'^ = 2a G'^ = a G* 24 = — 1 — 1 ^3 4 = ^ 


A 4 ,12 


C| 3 = 1 G| 3 = 1 G^ 4 = — 1 G 24 = 1 



Table 1: The structure constants of all 4-dim, real, Lie Algebras 
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5 Conclusion 



We have shown that the usage of the automorphism group is a very efficient way of 
identifying the true gravitational degrees of freedom for a simply transitive spatially 
homogeneous vacuum geometry. Many investigations have suffered from the failure of 
identifying these. In particular, if we wish to find the general solution under a given 
set of assumptions, then it is essential to ab initio identify the number of true degrees 
of freedom. At this point is we deem as appropriate to state that the Time-Dependent 
A.I.D.s were not only used to derive the second counting algorithm, but also to find some 
of the solutions exhibited in section 4. 

In this paper we specifically used this method to find the dimension of the set of all 
Ricci-flat spatially homogeneous models of dimension 4+1. Our main results are given 
in tables |2] and El 

Inspecting tables El and El it is seen that the most general types have 11 essential 
constants. Hence, in order to specify a certain solution under the above assumptions, 
we need to specify up to 11 parameters. The maximal number of parameters happens 
for the following two types 

A 3j8 © A 1 = VIII © R A 3t9 © A 1 = IX © R. 

Interestingly, these two algebras are the trivial extensions of the Bianchi type Lie algebras 
VIII and IX and not some indecomposable ones -as one might have expected. It is also 
noteworthy that, the set of the allowed numbers of the Essential Constants does not 
contain the numbers 1,3,4 and 5. This does not occur in 3+1 dimensions where the 
various models saturate all the range of values between 1 and 4. There, the models 
with the minimum number of essential constants are the Kasner (Type I) and Joseph 
(Type V). The corresponding 4+1 counterpart of Type I, i.e. 4Ai algebra is seen - 
by means of the algorithm- to contain 2 essential constants. Thus why the number 1 
is excluded. In fact this "hole" increases with the dimension, since the corresponding 
abelian types, will have d-2 essential constants, in d+1 dimensions. On the other hand, 
the 4+1 counterparts of the next "minimal" 3+1 models (Type V and its "neighbour" 
Type II with 1 and 2 essential constants respectively) i.e. the algebras A S 3 © A\ and 
© A\ have both 6 essential constants. The reason for this is that the number of 
the "would be constants" depend not only on the more components of the scale factor 
matrix 7 Q/ g (t) but also on the number of the linear constraints (the last being depended 
on the algebra). Thus from 2 the number of essential constants is lifted up to 6. Thus 
why the numbers between them i.e. 3,4,5 are also excluded. This sort of "irregularity" 
does not obtain for the rest of the cases, and thus all the numbers from 6 to 11 appear. 

We have also given some exact solutions, some of which are believed to be new. Only 
in two cases (4A\ and A 4 ' 5 ) the posited line element is the most general one. For the 
remaining types only special solutions are known. However, some of them - like the self- 
similar ones - may serve as asymptotes for more general solutions (this does, however, 
require a stability analysis within the class under consideration which to date is only 
done for the plane- wave solutions [20]). 
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TABLE 1 



Lie Algebra 




# of independent 
Linear Constraints 


dimAut(A) 


Essential 
Constants 


AA l 


20 





16 


2 


A 2 © 2Ai 


20 


4 


8 


6 




20 


4 


4 


10 


A 3 ,l © Ai 


20 


2 


10 


6 


A 3 , 2 © Ai 


20 


4 


6 


8 


A 3 ,3 © A x 


20 


4 


8 


6 


A 3 , 4 © A 1 


20 


3 


6 


9 


A% 5 ®A U < |a| < 1 


20 


4 


6 


8 


A 3 , 6 © A x 


20 


3 


6 


9 


A% 7 ®A 1 ,0<a 


20 


4 


6 


8 


A 3 , 8 © Ai and A 3<9 © 


20 


3 


4 


11 


A 4 ,i 


20 


3 


7 


8 


A£ 2 , a ^ {0, 1} 


20 


4 


6 


8 


A 4 , 2 


on 
20 


4 


Q 

o 





A 4 ,3 


on 
20 


4 


d 



o 
O 


A 4 ,4 


on 
20 


4 





o 
O 


A?|, a,/?e[-l,l)-{0},a^/J 


on 


4 





o 
O 


A 4 , 5 , « e 




on 
20 


4 


Q 

o 





A 4 ; 5 , a e 


-i,i) -{0} 


on 

20 


4 


8 





^4,5 


20 


4 


12 


2 


A%1 a^0,/?>0 


20 


4 


6 


8 


A 4 , 7 


20 


4 


5 


9 


A 4 ,s 


20 


3 


5 


10 


< 9 , < < 1 


20 


4 


5 


9 


4 1 

^4,9 


20 


4 


7 


7 


^4,9 


20 


4 


5 


9 


A 4 ,io 


20 


3 


5 


10 


A£ n , < a 


20 


4 


5 


9 


A 4 ,12 


20 


4 


4 


10 



Table 2: Essential Constants of 4+1 Spatially Homogeneous Models 
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TABLE 2 



Lie Algebra 




# of independent 
Linear Constraints 


dimAut(A) 


Essential 
Constants 


A% 5 © A u < |a| < 1 
' for a = -1/2, 


20 


3 


6 


9 


A% 2 , for a - -1, -3 


20 


3 


6 


9 


AH a ,/?e[-l,l)-{0} 
for 1 + 2a + /? = 
or 1 + 2/3 + a = 


20 


3 


6 


9 


^aeKi)-{0} 

for a = — 1 


20 


3 


8 


7 


^«e[-M)-{0} 

for ck = —1/3 


20 


2 


8 


8 


— 1 

<; 5 >e[-M)-{0} 

for a = —1 


20 


3 


8 


7 


A$>^0,/3>0 
for a = — (5 


20 


3 


6 


9 



Table 3: Essential Constants of 4+1 Spatially Homogeneous Exceptional Models 
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